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Measurements are presented from proton-proton collisions at centre-of-mass energies 
of ^/s = 0.9, 2.36 and 7 TeV recorded with the ATLAS detector at the LHC. Events were 
collected using a single-arm minimum-bias trigger The charged-particle multiplicity, its 
dependence on transverse momentum and pseudorapidity and the relationship between the 
mean transverse momentum and charged-particle multiplicity are measured. Measurements 
in different regions of phase-space are shown, providing diffraction-reduced measurements 
as well as more inclusive ones. The observed distributions are corrected to well-defined 
phase-space regions, using model-independent corrections. The results are compared to 
5^ I each other and to various Monte Carlo models, including a new AMBTl pythia6 tune. In 

all the kinematic regions considered, the particle multiplicities are higher than predicted by 
the Monte Carlo models. The central charged-particle multiplicity per event and unit of 
pseudorapidity, for tracks with pj > 100 MeV, is measured to be 3.483 ± 0.009 (stat) ± 
0.106 (syst) at ^ 0.9 TeV and 5.630 ± 0.003 (stat) ± 0.169 (syst) at = 7 TeV. 
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Additional Phase-Space Regions 



1 Introduction 



Inclusive charged-particle distributions have been previously measured in pp and pp collisions at a range 
of different centre-of-mass energies IllfTTll. These measurements provide insight into the strong interac- 
tions at low energy-scales. Several QCD-inspired models have been developed to interpret them. These 
models are frequently cast into Monte Carlo simulations with free parameters that can be constrained by 
measurements such as minimum bias distributions. These measurements contribute to the understand- 
ing of soft QCD; moreover, they are important to determination of biases on high-px phenomena due 
to underlying events and event pileup effects and are therefore of growing importance for future LHC 
physics. The measurements presented in this paper implement a similar strategy to that in [1]. A single- 
arm trigger overlapping with the acceptance of the tracking volume is used. Results are presented as 
inclusive-inelastic distributions, with minimal model-dependence; a minimum number of charged parti- 
cles within well-defined pj and 77 selection are required. 

This paper reports on measurements of primary charged-particle multiplicity distributions using the 
first ~190 ph^^ of data recorded by the ATLAS experiment at 7 TeV and ~1 p.b~^ at 0.9 TeV. At ^/s = 
0.9 TeV the sample is similar to that used for the first ATLAS minimum-bias publication Results 
are also presented at = 2.36 TeV where the track reconstruction setup differs significantly from that 
at the other energies, due to the Silicon Tracker (SCT) not being at nominal voltage. The integrated 
luminosity at this energy is estimated to be ~0.1 pb^^. 

The following distributions are measured in this paper: 

1 dA^ch 1 1 d^jVeh 1 dA^ev . , , 
A^ev d?7 A^ev 2npT drjdpT N^v dwch 

where pj is the charged particle momentum component transverse to the beam direction 0, 77 is the 
pseudorapidity of the particle, rich is the number of charged particles in an event, A'^ev is the number 
of events with a minimum number of charged particles within the selected kinematic range, A'^ch is the 
total number of charged particles in the data sample and (pj) is the average pj for a given number 
of charged particles q Primary charged particles are defined as charged particles with a mean lifetime 
T > 0.3 • 10~'° s either directly produced in pp interactions or from subsequent decays of particles with 
a shorter lifetime. 

The charged-particle multiplicity results are compared to particle level Monte Carlo (MC) predic- 
tions. Three different phase-space regions are considered in this paper, with varying selection both on 
the Pj and the number of charged particles per event; all phase-space regions require tracks within 
\r]\ < 2.5. Diffractive physics is expected to contribute mostly at low numbers of charged particles and 
at low track momentum. Therefore varying the selection on rich and pj in effect varies the relative con- 
tribution from diffractive events. Appendix IB] shows the results for two additional phase-space regions 
useful for Monte Carlo tuning. This measurement, with refined corrections and systematic uncertainty 
determination supersedes the results presented in [11 . 



2 The ATLAS Detector 

The ATLAS detector lITSl at the Large Hadron Collider (LHC) [19] covers almost the whole solid angle 
around the collision point with layers of tracking detectors, calorimeters and muon chambers. It has been 

'The ATLAS reference system is a Cartesian right-handed co-ordinate system, with the nominal collision point at the origin. 
The anti-clockwise beam direction defines the positive z-axis, while the positive j-axis is defined as pointing from the collision 
point to the centre of the LHC ring and the positive !/-axis points upwards. The azimuthal angle is measured around the beam 
axis and the polar angle 6 is measured with respect to the z-axis. The pseudorapidity is defined as = - lntan(£'/2). 

-The factor 2npj in the pj spectrum comes from the Lorentz invariant definition of the cross section in terms of d^p. Our 
results could thus be interpreted as the massless approximation to d^p. 
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designed to study a wide range of physics topics at LHC energies. For the measurements presented in 
this paper, the tracking devices and the trigger system are of particular importance. 

The ATLAS Inner Detector (ID) has full coverage in (p and covers the pseudorapidity range \r]\ < 2.5. 
It consists of a silicon pixel detector (Pixel), a silicon microstrip detector (SCT) and a transition radia- 
tion tracker (TRT). These detectors cover a sensitive radial distance from the interaction point of 50.5- 
150 mm, 299-560 mm and 563-1066 mm, respectively, and are immersed in a 2 T axial magnetic field. 
The inner-detector barrel (end-cap) parts consist of 3 (2x3) Pixel layers, 4 (2x9) double-layers of single- 
sided silicon microstrips with a 40 mrad stereo angle, and 73 (2x160) layers of TRT straws. Typical 
position resolutions are 10, 17 and 130 //m for the R-<p co-ordinate and, in case of the Pixel and SCT, 
115 and 580 /im for the second measured co-ordinate. A track from a charged particle traversing the 
barrel detector would typically have 1 1 silicon hits [1(3 pixel clusters and 8 strip clusters) and more than 
30 straw hits. 

For the runs at = 2.36 TeV, stable beams were not declared by the LHC; the high voltage on the 
SCT detector was thus not turned up to its nominal operating voltage but was left in standby mode. The 
Pixel detector was at nominal conditions for these runs. The hit efficiency in the SCT is thus significantly 
lower and special track reconstruction algorithms aie needed; the single hit efficiency at nominal voltage 
in the SCT barrel is above 99.7% EOl . while in standby it drops to ~ 60% for tracks perpendicular to the 
silicon surface. 

The ATLAS detector has a three-level trigger system: Level 1 (LI), Level 2 (L2) and Event Fil- 
ter (EF). For this measurement, the trigger relies on the LI signals from the Beam Pickup Timing de- 
vices (BPTX) and the Minimum Bias Trigger Scintillators (MBTS). The BPTX stations are composed 
of electrostatic button pick-up detectors attached to the beam pipe at ±175 m from the centre of the 
ATLAS detector. The coincidence of the BPTX signal between the two sides of the detector is used to 
determine when bunches are colliding in the centre of the ATLAS detector. The MBTS are mounted at 
each end of the detector in front of the liquid-argon end-cap calorimeter cryostats &tz - +3.56 m. They 
are segmented into eight sectors in azimuth and two rings in pseudorapidity (2.09 < \q\ < 2.82 and 
2.82 < \ri\ < 3.84). Data were collected for this analysis using a trigger requiring a BPTX coincidence 
and MBTS trigger signals. The MBTS trigger used for this paper is configured to require one hit above 
threshold from either side of the detector, referred to as a single-arm trigger. The efficiency of this trig- 
ger is studied with a separate prescaled LI BPTX trigger, filtered to obtain inelastic interactions by Inner 
Detector requirements at L2 and EF, the latter only for the 900 GeV data. 

3 Monte Carlo Simulation 

Inclusive minimum bias data are modelled using three components in the pythia6 ll2Tll Monte Carlo 
(MC) event generator: non-diffractive (ND), single- (SD) and double-diffractive (DD). Non-difFractive 
processes are modelled from two-to-two processes as described in this section. Diffractive process mod- 
elling is described in Sec. 13. II 

Low-/7t scattering processes may be described by lowest-order perturbative Quantum Chromody- 
namics (QCD) two-to-two parton scatters, where the divergence of the cross section at /jt = is regulated 
by phenomenological models. The pythia6 MC event generator implements several of these models. The 
parameters of these models have been tuned to describe charged-hadron production and the underlying 
event in pp and pp data at centre-of-mass energies between 200 GeV and 1.96 TeV. 

Samples of MC events were produced for single-difFractive, double-diffractive and non-diffractive 
processes using the pythia6 generator 0. The ATLAS MC09 pythia tune f2T\ uses a specific set of 

'a hit is a measurement point assigned to a track. 
PYTHIA version 6.4.21 
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optimised parameters; it employs the MRST LO* parton density functions (PDFs) and the pj- 
ordered parton shower |24|. A tune is a particular configuration or set of values of the parameters of 
the particular Monte Carlo model. These parameters were derived by tuning to the underlying event 
(UE) and minimum-bias data from the Tevatron at 630 GeV to 1.96 TeV. The MC samples generated 
with this tune are used to determine detector acceptances and efficiencies and to correct the data. MC 
samples were produced at all three centre-of-mass energies considered in this paper. The non-diffractive, 
single-difFractive and double-diffractive contributions in the generated samples are mixed according to 
the generator cross sections. 

All the events are processed through the ATLAS detector simulation program Il25l . which is based 
on geant4 1261 . They are then reconstructed and analysed by the same program chain used for the 
data. Particular attention was devoted to the description in the simulation of the size and position of 
the collision beam spot and of the detailed detector conditions during data taking. The MC09 pythia6 
samples are used to derive the detector corrections for these measurements. The MC samples at 2.36 TeV 
were generated assuming nominal detector conditions. 

For the purpose of comparing the present measurements to different phenomenological models de- 
scribing minimum-bias events, the following additional particle level MC samples were generated: 

• the new ATLAS Minimum Bias Tune 1 (AMBTl) pythia6 tune described in Sec. 13.21 

• the DW 1 27 1 pythia6 tune, which uses virtuality-ordered showers and was derived to describe the 
CDF Run II underlying event and Drell-Yan data; 

• the pythiaS generator I EH, in which the diffraction model produces much harder pj and rich 
spectra for the single- and double-diffractive contributions than pythia6. The default parton shower 
model is similar to the one used in pythia6 MC09; 

• the PHOJET generator H [29], which is used as an alternative model to PYiHiA-based generators. 
PHOJET relies on PYTHiA6L|for the fragmentation of partons. 

3.1 Diffractive Models 

PYTHIA6, pythiaS and phojet model the diffractive components very differently. Here we mostly describe 
the model implemented in pythia6. The pythia6 diffraction is based on a Regge-based pomeron model 
to generate the cross-section and generate the diffractive mass and momentum transfer 113011311 . To allow 
the Regge model to cover the full phase-space, empirical corrections are introduced fTP\. These have 
the effect of enhancing the production of small masses and suppressing production near the kinematic 
limit. Particle production from low mass states (Mx < 1 GeV) is treated as an isotropic two body decay. 
Particle production from high mass states is based on the string model. Two string configurations are 
possible depending on whether the pomeron couples to a quark or gluon |[2T1l . 

The pythiaS model uses the same model as pythia6 to generate the cross-section and generate the 
diffractive mass and momentum transfer. The particle production for low mass states uses the string 
model but for higher masses {Mx > 10 GeV) a perturbative element based on pomeron-proton scattering 
is introduced. The non-perturbative string model introduces a mass dependence on the relative probabil- 
ity of the pomeron scattering off a quark to scattering off a gluon, which enhances the gluon probability 
at high masses. The perturbative pomeron-proton scattering uses HERA diffractive PDFs |[32l and the 
standard multiple interactions framework is used to generate the parton-parton scattering. The introduc- 
tion of the perturbative pomeron-proton scattering results in a harder pj and multiplicity spectrum for 

'pYTHiA version 8.130 
*PHOJET version 1.12.1.35 
'pYTHiA version 6.1.15 
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diffractive events generated with pythiaS compared to those generated with pythia6 However, it 
should be noted that relatively little tuning has been made of the diffractive processes in pythia6 and 
pythiaS. 

PHOJET is based on the dual parton model. It generates a harder pj and multiplicity spectrum in 
diffractive events than pythia6. The new diffraction model of pythiaS generates distributions quite similar 
to those from phojet |[33l . 

3.2 PYTHIA 6 ATLAS Minimum Bias Tune 1 

Before the start of the LHC, an ATLAS tune to pythia6 with MRST LO* PDFs using Tevatron un- 
derlying event and minimum bias data was produced, the so-called MC09 tune [22|. The first ATLAS 
measurements of charged particle production at the LHC [1] measured the charged particle production 
at = 0.9 TeV in the central region to be 5-15% higher than the Monte Carlo models predict. In addi- 
tion, neither the high rich nor the high pj distributions were well described by this tune and the (pj) was 
overestimated in events with rich > 20. A new tune, AMBTl, was developed in order to adapt the free 
parameters of the non-diffractive models to the new experimental data at = 0.9 TeV and ^/s = 7 TeV, 
using the same PDFs and pythia6 model choices as MC09. 

The AMBTl tune is obtained by tuning to ATLAS minimum bias data at both ^fs = 0.9 TeV and 
^/s = 7 TeV in a diffraction-reduced phase-space that is presented in this paper: rich > 6, pi > 500 MeV, 
I77I < 2.5. The tune was derived using preliminary versions of these distributions [34]. The starting point 
for this tune is the ATLAS MC09c 1221 pythia6 tune. MC09c is an extension of the ATLAS MC09 tune 
where the strength of the colour reconnection (CR) was tuned to describe the (pj) vs. rich distributions 
measured by CDF in pp collisions at the Tevatron [Tl. 

Charged particle distributions are sensitive to multi-parton interactions (MPI) and colour reconnec- 
tion of the hadronic final state |[35l : the MPI are regulated by a low pj cut-off and the matter overlap 
distribution of the two protons in which the additional partonic scattering takes place. These are the 
main parameters varied for this new tune. Parameters related to final state radiation, hadronisation and 
fragmentation are not tuned, as these are constrained by many LEP results. No changes to the diffraction 
model are made. The model parameters are adapted in order to best describe these new distributions over 
the full range while maintaining consistency with the Tevatron results. For the data MC comparisons 
the Rivet 1^ |36| package is used; the tuning is done using the professor package 

I^ISIIMI- Tabled] 

summarizes the parameters varied in this tune; the meaning of the parameters are given below. 

MPI Parameters The size of the MPI component in the pythia6 model is regulated by a simple cut-off 
parameter for the pj of two-to-two scattering processes. This cut-off parameter is fixed at a reference 
energy, which is generally taken as 1.8 TeV. The cut-off at this reference scale is called PARP(82). It 
is then rescaled for other centre-of-mass energies using a parameter PARP(90). The rescaling is done 
according to the following formula: 



The amount of scattering is described by the matter overlap distribution between the two protons, 
which regulates how many central, hard scatterings and how many less central, softer scatterings occur. 
This distribution is modelled as a double Gaussian probability density function. The parameter PARP(83) 
describes the fraction of matter in the naiTower of the two Gaussian functions. The size of this narrower 
Gaussian is given as a fraction PARP(84) of the wider, main radius. The optimal value for this parameter 

version 1.2.2a0 
'version l.O.OaO 




PARP(90) 



(1) 
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was found in a first tuning run. Further variations of the matter fraction in the narrower cone were found 
to not have a significant influence on the main distributions used for tuning. 

Colour Reconnection Parameters The colour reconnection scenario of pythia used in MC09c min- 
imises the total string length between partons. The probability that a given string piece does not partici- 
pate in the CR is given by (1 - PARP(78))"^' , where nui is the number of multi-parton interactions 1211 : 
the larger the parameter, the smaller the probability of the string piece not participating. In addition to this 
parameter, an additional parameter PARP(77) is present in pythia; it is used to describe a suppression fac- 
tor for the CR of fast moving string pieces. The suppression factor is given by 1/(1 -i- PARP(77)^ • Pavg)> 
where pl^^ is a measure of the average squared momentum that hadrons produced by the string piece 
would have. 

Additional Parameters Investigated In an initial study, the cut-off parameter for initial state radiation 
(PARP(62)) and the cut-off for momentum smearing in primordial kj_ (PARP(93)) were considered. The 
optimal values for these parameters were found in a first tuning run, further variation of those parameters 
was not found to have a significant influence on the main distributions used for tuning. 

Distributions Used The tune described in this paper focuses on the ATLAS minimum bias data. It 
primarily attempts to improve the description of the high pj and high n^h distributions observed. For 
the pt spectrum, only particles above 5 GeV are considered. For the rich spectrum, only events with 
20 or more tracks are used in the tune. For the (pj) vs. rich distribution, only events with ten or more 
tracks are considered. The full r] distribution is used. For completeness, the preliminary underlying event 
results II39II40II are included in the plateau region; however, due to the limited statistics, these data have 
only very small impact on the tune. 

Tevatron data in the energy range of 630 GeV to 1 .96 TeV are included in the tune, but with a weight 
which is ten times lower than that of the ATLAS data. This weighting allows a check of the consistency 
of the resulting tune with the Tevatron data while forcing the ATLAS data to drive the tuning process. 
Similar datasets were used for the MC09c tune. The charged particle multiplicity shown in ll4T1l was not 
included in the tune as no variation of the tuning parameters considered was able to fit both the ATLAS 
and the CDF distributions simultaneously. App. lAl shows a full list of the distributions and the ranges 
considered by the tune. 



Results The final parameter values resulting from the tune are shown in Table [T] 



Parameter 


Related model 


MC09c value 


scanning range 


AMBTl value 


PARP(90) 


MPI (energy extrapolation) 


0.2487 


0.18-0.28 


0.250 


PARP(82) 


MPI (/J™") 


2.31 


2.1-2.5 


2.292 


PARP(84) 


MPI matter overlap (core size) 


0.7 


0.0-1.0 


0.651 


PARP(83) 


MPI matter overlap (fraction in core) 


0.8 


fixed 


0.356 


PARP(78) 


CR strength 


0.224 


0.2 - 0.6 


0.538 


PARP(77) 


CR suppression 


0.0 


0.25-1.15 


1.016 


PARP(93) 


Primordial kj_ 


5.0 


fixed 


10.0 


PARP(62) 


ISR cut-off 


1.0 


fixed 


1.025 



Table 1: Comparison of MC09c and AMBTl parameters. The ranges of the parameter variations scanned 
are also given. The parameters declared as 'fixed' were fixed to the values obtained after an initial pass 
of the tuning. 
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4 Data Selection 



Events in which the Inner Detector was fully operational and the solenoid magnet was on are used for 
this analysis for both ^/s = 0.9 TeV and ^/s = 1 TeV. During this data-taking period, more than 97% 
of the Pixel detector, 99% of the SCT and 98% of the TRT were operational. At = 2.36 TeV the 
requirements are the same, except for the SCT being in standby. 

Events were selected from colliding proton bunches in which the MBTS trigger recorded one or 
more counters above threshold on either side. The maximum instantaneous luminosity is approximately 
1.9 X 10^^ cm~^ s~' at 7 TeV. The probability of additional interactions in the same bunch crossing is 
estimated to be of the order of 0.1%. In order to perform an inclusive-inelastic measurement, no further 
requirements beyond the MBTS trigger are applied. 

In order to better understand the track reconstruction performance at ^/s = 2.36 TeV, during which 
time the SCT was in standby, additional data at = 0.9 TeV were taken with the SCT in standby for 
part of a run. This enables the derivation of data-driven corrections to the track reconstruction efficiency, 
as described in Sec. 16.41 

4.1 Different Phase-Space Regions Considered 

Three separate phase-space regions are considered in the main part of this paper with varying contribu- 
tions from diffractive events: 

• at least one charged particle in the kinematic range I77I < 2.5 and pj > 500 MeV, 

• at least two charged particles in the kinematic range \t]\ < 2.5 and /jt > 100 MeV, 

• at least six charged particles in the kinematic range \t]\ < 2.5 and pj > 500 MeV. 

The first of these phase-space regions is studied at all three centre-of-mass energies. This is the 
region that allows us to best investigate the evolution of charged-multiplicity distributions as a function 
of centre-of-mass energy and thus constrain the MC parameters that dictate the energy extrapolation of 
the models. The second measures the most inclusive charged-particle spectra and is also used as the 
basis for the model-dependent extrapolation to pj = 0; in this phase-space region results at = 0.9 
and 7 TeV are shown. The third phase-space region considered is similar to the first but with a higher cut 
on the number of charged particles, thus reducing the expected contribution from diffractive events in the 
sample. These distributions are measured for both 0.9 and 7 TeV. This is the phase-space region which 
was used to produce the new AMBTl tune. At 2.36 TeV only the first phase-space region is measured. 
Two additional phase-space regions are presented in App.lBl 

The relative contribution from diffractive events varies widely between Monte Carlo models and 
depends strongly on the phase-space region selection applied. The diffractive contribution is constrained 
very little by previous data. Table |2] shows the predicted fractions of simulated events originating from 
diffractive processes, as predicted by pythia6, pythiaS and phojet; the values for the different tunes of 
pythia6 are found to be similar because the acceptances of the different non-diffractive models do not 
change significantly and the diffractive models are identical. The large difference in predictions between 
the models is one of the motivations for not making any model-dependent corrections to the experimental 
data, as such corrections would vary significantly depending on which MC model is used to derive them. 

4.2 Event Selection 

To reduce the contribution from background events and non-primary tracks, as well as to minimise the 
systematic uncertainties, the events are required to satisfy the following criteria: 
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Phase-Space Region 


= 0.9 TeV 


= 7 TeV 




min «ch 


min pj (MeV) 


pythia6 


pythiaS 


PHOJET 


pythia6 


pythiaS 


PHOJET 


2 


100 


22% 


22% 


20% 


21% 


21% 


14% 


1 


500 


16% 


21% 


19% 


17% 


21% 


14% 


6 


500 


0.4% 


5% 


8% 


0.4% 


10% 


8% 



Table 2: Fraction of simulated events originating from diffractive processes, as predicted by pythia6, 
pythia8 and phojet in the three phase-space regions measured in this paper at both ^fs = 0.9 TeV and 
^/s = 1 TeV. All results are for I77I < 2.5. 



• to have triggered the single-arm, single-counter level 1 minimum bias trigger scintillators 

• the presence of a primary vertex 14211 reconstructed using the beam spot information H3l and at 
least two tracks, each with: 

- Pt> 100 MeV, 

- a transverse distance of closest approach with respect to the beam-spot position Id^^ \ < 4 mm; 

• the rejection of events with a second vertex containing four or more tracks, to remove events with 
more than one interaction per bunch crossing; 

• a minimum number of tracks, depending on the particular phase-space region, as described in 
Sec. 1431 

4.3 Track Reconstruction Algorithms 

Tracks are reconstructed offline within the full acceptance range |?7| < 2.5 of the Inner Detector f44','45l. 
Track candidates are reconstructed by requiring a minimum number of silicon hits and then extrapolated 
to include measurements in the TRT. Due to the SCT being in standby mode at 2.36 TeV, different track 
reconstruction algorithms are needed; at 0.9 and 7 TeV, the reconstruction algorithms are collectively 
referred to as full tracks. The analysis at = 2.36 TeV has been performed using two complementary 
methods for reconstructing tracks. The first reconstructs tracks using pixel detector information only, 
denoted Pixel tracks. The second uses tracks reconstructed from the full Inner Detector information, 
denoted ID tracks 0. 

4.3.1 Algorithms for 0.9 and 7 TeV 

For the measurements at 0.9 and 7 TeV, two different track reconstruction algorithms are used. The 
algorithm used for the previous minimum-bias publication [ 1 j is used with a lower pj threshold cut at 
100 MeV. An additional algorithm configuration is run using only the hits that have not been used by 
the first algorithm. This additional algorithm uses wider initial roads and has a looser requirement on 
the number of silicon hits. This second algorithm contributes around 60% of the tracks from 100 to 
150 MeV, mostly due to the tracks having too low a momentum to go far enough in the SCT detector to 
satisfy the silicon hit requirement of the original algorithm; this fraction decreases rapidly, reaching less 
than 2% at 200 MeV. 

Tracks are required to pass the selection criteria shown in Table [3j the column labelled Full Tracks 
refers to the algorithms used at 0.9 and 7 TeV. The transverse, do, and longitudinal, z.o, impact parameters 

'In the context of the other analyses, ID tracks are referred to as track for brevity. 
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a/v — Q and 7 TeV 


= 2.36 TeV 




Pull Tr^ipV^ 


ID Tracks 


Pixel Tracks 


> 1 on or 500 MeV 


YES 


YES 


YES 


< 2.5 


YES 


YES 


YES 


layer-0 hit if expected 


YES 


YES 


YES(*) 


> 1 Pixel hit 


YES 


YES 


YES 


> 2, 4 or 6 SCT hits for tracks (**) 


YES 


NO 


NO 


l^ol < 1-5 mm and |zol • sin0 < 1.5 mm 


YES 


YES 


YES(***) 


probability > 0.01 for pj > 10 GeV 


YES 


N/A 


N/A 



Table 3: Selection criteria applied to tracks for the full reconstruction, ID tracks and Pixel tracks. The 
transverse momentum cut applied depends on the phase-space region in question. (*) For the Pixel track 
method the layer-0 is required even if not expected. (**) The SCT hit selection are for p-Y < 200, 
200 < Pi < 300 or pi > 300 MeV, respectively. (***) For the Pixel track method, the and zq 
selection are after the track refitting is performed (see Sec. 14.3.21 . 



are calculated with respect to the event primary vertex. The layer-0 selection requires a hit in the inner- 
most layer of the Pixel detector if a hit is expected H. The track-fit probability H cut is applied to 
remove tracks with mis-measured pi due to mis-alignment or nuclear interactions. 

These tracks are used to produce the corrected distributions and will be referred to as selected tracks. 
The multiplicity of selected tracks within an event is denoted by «sei- The tracks used by the vertex 
reconstruction algorithm are very similar to those used for the analysis; the p^ threshold is also 100 MeV. 
Due to the requirement that the vertex be made from a minimum of two such tracks and the fact that we do 
not wish to correct our measurement outside of the observed phase-space region, the minimum number 
of particles per event for the phase-space region with px > 100 MeV also needs to be set at two. Table|4] 
shows the total number of selected events and tracks for all phase-space regions considered. 

Trigger and vertex reconstruction efficiencies are parameterised as a function of n^^^. n^^^ is defined 
as the number of tracks passing all of the track selection requirements except for the constraints with 
respect to the primary vertex; instead, the unsigned transverse impact parameter with respect to the beam 
spot, \d^^\, is required to be less than 1.8 mm. 



Phase-Space Region 


Vi = 


0.9 TeV 




7 TeV 




2.36 TeV 




min pi 


Full Tracks 


Full Tracks 


ID Tracks (Pixel Tracks) 




(MeV) 


Events 


Tracks 


Events 


Tracks 


Events 


Tracks 


2 


100 


357,523 


4,532,663 


10,066,072 


209,809,430 






1 


500 


334,411 


1,854,930 


9,619,049 


97,224,268 


5,929 (5,983) 


38,983 (44,788) 


6 


500 


124,782 


1,287,898 


5,395,381 


85,587,104 







Table 4: Number of events and tracks in the three phase-space regions at each centre-of-mass energy 
considered in this paper. 



"A hit is expected if the extrapolated track crosses an active region of a Pixel module that has not been disabled. 

^^This probability function is computed as 1 - P(nii„f/2,x^/2), where P(ni„f/2,x^/2) is the incomplete gamma function and 
«dof is the number of degrees of freedom of the fit. It represents the probability that an observed exceeds the observed value 
for a correct model. 
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4.3.2 Track Reconstruction Algorithms at 2.36 TeV 

Operation of the SCT at standby voltage during 2.36 TeV data taking led to reduced SCT hit efficiency. 
Consequently, ID tracks are reconstructed at this centre-of-mass energy using looser requirements on 
the numbers of hits and holes Il44ll45l . There are no simulation samples that fully describe the SCT 
operating at reduced voltage. A technique to emulate the impact of operating the SCT in standby was 
developed in simulation; this corrects the Monte Carlo without re-simulation by modifying the silicon 
clusterisation algorithm used to study the tracking performance. However, the final ID track efficiency 
at = 2.36 TeV was determined using a correction to the track reconstruction efficiency derived from 
data at = 0.9 TeV. 

Pixel tracks were reconstructed using the standard track reconstruction algorithms limited to Pixel 
hits and with different track requirements. There is little redundant information, because at least three 
measurement points are needed to obtain a momentum measurement and the average number of Pixel 
hits per track is three in the barrel. Therefore the Pixel track reconstruction efficiency is very sensitive 
to the location of inactive Pixel modules. The total distance between the first and the last measurement 
point in the pixel detector, as well as the limited number of measurement points per track, limit the 
momentum resolution of the tracks; therefore the Pixel tracks were refit using the reconstructed primary 
vertex as an additional measurement point. The refitting improves the momentum resolution by almost 
a factor of two. However, the Pixel track momentum resolution remains a factor of three worse than the 
resolution of ID tracks. 

The selection criteria used to define good Pixel and ID tracks are shown in Table [3l The total number 
of accepted events and tracks at this energy are shown in TablelH These two track reconstruction methods 
have different limitations; the method with the best possible measurement for a given variable is chosen 
when producing the final plots. The Pixel track method is used for the rich and rj distributions, while the 
ID track method is used for the pj spectrum measurement; the (pj) distribution is not produced for this 
energy as neither method is able to describe both the number of particles and their pj accurately. 

5 Background Contribution 

5.1 Event Backgrounds 

There are three possible sources of background events that can contaminate the selected sample: cosmic 
rays, beam-induced background and the presence of another collision inside the same bunch crossing. 
The fraction of cosmic ray background events was estimated in |1], where it was found to be smaller 
than 10"^. Beam-induced backgrounds are estimated from non-colliding empty bunches using the same 
method as described in IT] ; after final event selection, fewer than 0. 1 % of events are predicted to originate 
from beam-induced backgrounds. The reconstructed primary vertex requirement is particularly useful in 
suppressing the beam-induced background. The instantaneous luminosity at ^/s = 7 TeV is high enough 
that the effect of multiple collisions inside the same bunch crossing cannot be ignored. Events are rejected 
if they have a second vertex with four or more tracks H. After this cut, the fraction of events with more 
than one interaction in the same bunch crossing is measured to be about 0.1%; the residual effect is thus 
neglected. At the lower centre-of-mass energies, the rate of multiple interactions is lower and thus also 
neglected. 

*A hole is defined as an absence of a hit when it is expected given the track trajectory. 

'Events with two vertices with fewer than four tracks are dominated by events where a secondary interaction is reconstructed 
as another primary vertex and are thus not removed from our data samples. 
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5.2 Backgrounds to Primary Tracks 



Primary charged-particle multiplicities are measured from selected-track distributions after correcting 
for the fraction of non-primary particles in the sample. Non-primary tracks are mostly due to hadronic 
interactions, photon conversions and decays of long-lived particles, as well as a small fraction of fake 
tracks. Their contribution is estimated using MC predictions for the shape of the do distribution for 
primaries, non-primaries from electrons and other non-primaries. The separation between non-primaries 
from electrons and non-electrons is needed as the electrons are mostly from conversions in the detector 
material and would thus be sensitive to a mis-modeling of the detector material, whereas the non-electron 
non-primary tracks are mostly from long-lived particles and this fraction is thus also sensitive to the 
underlying physics. The Gaussian peak of the do distribution, shown in Fig. [T]for 100 < pj < 150 GeV, 
is dominated by the primary tracks and their resolution. The non-primary tracks populate the tails. The 
dominant contribution to non-primary tracks inside the acceptance cut on |<io| comes from non-electrons. 

The primary, electron non-primary and non-electron non-primary do distributions are obtained from 
MC and used as templates to extract the relative fractions in data. A fit is performed in the side-bands of 
the distribution, i.e. outside the range in do used for selecting tracks. The fractions of primary, electron 
non-primary and non-electron non-primary tracks are all allowed to float with the total number of events 
constrained to that of the data. The contribution of non-primaries from electrons within the analysis 
acceptance of 1.5 mm is small, while it dominates at high values of \do\. The requirement on having a 
hit on layer-0 suppresses this contribution enough to allow the fit to be performed down to the lowest 
Pi region. The fit is performed in bins of 50 MeV in pj from 100 to 500 MeV. A single fit is used 
for all tracks with pj > 500 MeV; in this bin the distinction is not made between the two sources of 
non-primary tracks. The fraction of non-primary tracks varies from 3.4% for 100 < pj < 150 MeV to 
1.6% above 500 MeV at ^/s = 7 TeV. Figure[T]shows the observed do distribution for the bin 100 < pj < 
150 MeV compared to the MC predictions after the fit. 

E 
E 

CJ 

o 

CO fi 

-g 10^ 



10^ 
10* 



-10 -8 -6 -4 -2 2 4 6 8 10 

do [n^fTi] 

Figure 1: Transverse impact parameter, do, distribution at ^/s = 7 TeV for primary (blue short dashed) 
and non-primary particles after scaling them to the best fit value for 100 < /jt < 150 MeV. The non- 
primary particles are split into electrons (pink long-dashed) and non-electrons (green dot-dashed). The 
full red curve shows the non-diffractive (ND) MC prediction for the sum over the three components 
which agrees well with the data (black points). 

Systematic Uncertainties The full difference between the non-primary fraction in MC and that in data 
obtained using the fit is taken as a systematic uncertainty. The largest difli'erence is found to be an increase 
of non-primaries in data by 25% relative to the MC for pj > 500 MeV. This conservative estimate is 
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taken to be constant as a function of pj and results in only a small effect, up to 0.9%, on the final 
corrected distributions. In order to estimate the effect of the choice of the variable used to obtain the fit, 
the fraction of primary and non-primary track contributions are obtained by fitting the zo distributions. 
The difference is measured to be 12% in the first bin, 8% in the last bin and less than 4% in all other 
bins; this difference is taken as a source of systematic uncertainty. The estimated number of non-primary 
tracks in \do\ < 1.5 mm is found to be stable with respect to a change in the fit range of 1 mm in all pj 
bins except the first one (100 < pj < 150 MeV), where a 10% difference is observed; this difference is 
taken as a systematic uncertainty. The fraction of non-primary tracks is found to be independent of ?isei> 
but shows a small dependence on 77, taken as a small systematic uncertainty of 0. 1%. 

The total uncertainty on the fraction of non-primary tracks is taken as the sum in quadrature of 
all these effects. The total relative uncertainty on the measured distributions at = 0.9 TeV and 
^/s ^ 7 TeV is 1.0% for the first pj bin, decreasing to 0.5% above 500 MeV. At = 2.36 TeV this 
uncertainty for the Pixel track method is 0.6%. 

6 Selection Efficiency 

The data are corrected to obtain inclusive spectra for charged primary particles satisfying the different 
phase-space region requirements. These corrections include inefficiencies due to trigger selection, vertex 
and track reconstruction. They also account for effects due to the momentum scale and resolution and 
for the residual background from non-primary tracks. 

In the following sections the methods used to obtain these efficiencies, as well as the systematic 
uncertainties associated with them are described. Plots are shown for the phase-space region rich ^ 2, 
Pt > 100 MeV, \r]\ < 2.5 at ^/s = 7 TeV, but similar conclusions can be drawn at the other energies and 
phase-space regions. 

6.1 Trigger Efficiency 

The trigger efficiency, etrig> is measured from a data sample selected using a control trigger. The control 
trigger used for this analysis selects events from random filled bunch crossings which are then filtered 
at L2. At = 0.9 TeV the L2 filter requires a minimum of seven pixel clusters and seven SCT hits 
and the EF requires at least one track with pj > 200 MeV. At ^fs = 7 TeV the L2 requirement is 
loosened to four pixel clusters and four SCT hits. No EF requirements are made at this energy. The 
vertex requirement for selected tracks is removed for these trigger studies, to account for correlations 
between the trigger and vertex reconstruction efficiencies. The trigger efficiency is determined by taking 
the ratio of events from the control trigger in which the LI MBTS also accepted the event, over the total 
number of events in the control sample. For ^/s = 2.36 TeV there is not sufficient data to measure the 
trigger efficiency and thus the = 0.9 TeV parametrisation is used to correct the 2.36 TeV data. 

The trigger efficiency is parametrised as a function of n^^^; it is 97% (99%) in the first n^^^ bin and 
rapidly increases to nearly 100% for nf^^ > 2, pj > 100 MeV (n^J > 1, Pt > 500 MeV). The trigger 
requirement is found to introduce no observable bias in the pj and rj distributions of selected tracks 
within the statistical uncertainties of the the data recorded with the control trigger. The resulting trigger 
efficiency is shown in Fig. |2^ for the phase-space region with >2,pT > 100 MeV at = 7 TeV. 

Systematic Uncertainties Since there is no vertex requirement in the data sample used to measure the 
trigger efficiency, it is not possible to make the same impact-parameter selection as is made on the final 
selected tracks. In order to study potential effects due to this, the trigger efficiency is measured after 
applying the impact-parameter constraints with respect to the primary vertex if available or with respect 
to the beam spot if not. The difference in the efficiency obtained this way and in the nominal way is 
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Figure 2: Trigger efficiency (a) and vertex reconstruction efficiency (b) with respect to the event selection, 



as a function of the number of reconstructed tracks before the vertex requirement (n^?). The track 



reconstruction efficiency as a function of rj (c) and pj (d) is derived from non-diffractive (ND) MC. The 
statistical errors are shown as black lines, the total errors as green shaded areas. All distributions are 
shown at ^/s = 1 TeV for rich > 2, pj > 100 MeV, I77I < 2.5. For the vertex and trigger efficiencies, the 
selection requires n^^^ > 2. 
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considered as a systematic uncertainty. This variation provides a conservative estimate of the effect of 
beam-induced background and non-primary tracks on the trigger efficiency at low values of n^^^. The 
systematic uncertainty arising from possible correlation of the MBTS trigger with the control trigger is 
studied using simulation, and the effect of correlations on the trigger efficiency is found to be less than 
0.1%. The total systematic uncertainty on the trigger efficiency determination, which also includes the 
statistical uncertainty on the control sample, is of the order of 1% in first bin, decreasing rapidly as 
n^^ increases. 

sel 

6.2 Vertex Reconstruction Efficiency 

The vertex reconstruction efficiency, £vtx, is determined from data by taking the ratio of triggered events 
with a reconstructed vertex to the total number of triggered events, after removing the expected contri- 
bution from beam background events. The efficiency is measured to be 90-92% in the first rP^^ bin for 
the different energies and phase-space regions; it rapidly rises to 100% at higher track multiplicities. The 
vertex reconstruction efficiency at ^/s = 1 TeV for n^^^ > 2, pj > 100 MeV is shown in Fig. [2j) as a 
function of n^? . 

sel 

The dependence of the vertex reconstruction efficiency on the rj and pj of the selected tracks is 
studied as well as the dependence on the projection along the beam-axis of the separation between the 
perigees |^f the tracks (Az), for events with more than one track. For all phase-space regions, only the 
dominant effect is corrected for as the other effect is always found to be significantly smaller and would 
thus not affect the final result. 

For the lower pj threshold selection, a strong dependence is observed as a function of Az for events 
with two tracks; this bias is corrected for in the analysis using two different parametrisations depending 
on the Pj of the lowest pj track: one for tracks below 200 MeV and one for those above that threshold. 
The dependence on the vertex reconstruction efficiency due to the 77 of the tracks is found to be smaller 
than the Az correction and is neglected for this phase-space region. For the 500 MeV pj threshold 
selection, the rj dependence is corrected for events with rfi^^ = 1 . For events with higher multiplicities 
the Az dependence is found to be very small and is neglected. 

Systematic Uncertainties The difference between the vertex reconstruction efficiency measured with 
beam background removal and the vertex reconstruction efficiency measured without beam background 
removal is assigned as the systematic uncertainty on the vertex reconstruction efficiency. For determi- 
nation of this difference, the contribution of beam-related backgrounds is estimated using non-colliding 
bunches, as in [ 1]. The highest rate of beam-related background is found in the phase-space region with 
Pt > 100 MeV at 900 GeV, where it is 0.8% without vertex selection and 0.2% with vertex selec- 
tion, although it is found to decrease rapidly at higher multiplicities. (This beam-related background 
contribution is larger than that given if Sec. |5] where a reconstructed primary vertex was required.) The 
total uncertainty due to the vertex reconstruction efficiency is significantly below 1 % for all phase-space 
regions at all energies. Fig|2j5 shows the total error for the phase-space region with pj > 100 MeV at 
V* = 7 TeV. 

6.3 Track Reconstruction Efficiency for the 0.9 and 7 TeV Data Samples 

The track reconstruction efficiency, Strk, determined from MC, is parametrised in bins of pj and rj. 
The excellent agreement between data and MC of basic track quantities for tracks above 500 MeV was 
previously demonstrated |T|. Figure[3]highlights the agreement for tracks in the additional range covered 
in this paper, \00 < pr < 500 MeV. 

*The perigee of a track is liere tlie point of closest approach of the track and the coordinate origin (0,0,0). 
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Figure 3: Comparison between data and simulation at ^/s = 1 TeV for tracks with transverse momentum 
between 100 and 500 MeV: the average number of silicon hits on reconstructed track as a function of 
77 in the SCT (a) and Pixel (b) detectors, the transverse impact parameter (c) and longitudinal impact 
parameter multiplied by sin 9 (d). The inserts for the impact parameter plots show the log-scale plots. 
The pj distribution of the tracks in non-diffractive (ND) MC is re-weighted to match the data and the 
number of events is scaled to the data. 



The track reconstruction efficiency is defined as: 

A^c"(;'T,'7) 



where pj and rj are generated particle properties, N^^^^^^{pt,t]) is the number of reconstructed tracks 
matched to a generated charged particle and Nge„{pT,T]) is the number of generated charged particles 
in that bin. The matching between a generated particle and a reconstructed track uses a cone-matching 
algorithm in the rj-cp plane, associating the particle to the track with the smallest AR = yJ{A(p)^ + {Ar])^ 
within a cone of radius 0.15. In addition, the particle trajectory must be compatible with the position of 
one of the pixel hits of the track. The larger cone size than in m is needed to account for the degraded 
resolution at lower track pj. 

The resulting reconstruction efficiency as a function of t] integrated over pj is shown in Fig. |2}: at 
yfs = 7 TeV for the phase-space region with the lowest pj threshold. The track reconsti"uction efficiency 
is lower in the region \r]\ > 1 due to particles passing through more material in that region. Figure |2}i 
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shows the efficiency as a function of pj integrated over 77. The initial rise with pi is due to the requirement 
on the minimum number of silicon hits required in the analysis, which indirectly constrains the tracks to 
pass through a minimum number of detector layers and thus have a minimum /jt- 

Systematic Uncertainties As the track reconstruction efficiency is determined from MC, the main 
systematic uncertainties result from the level of agreement between data and MC. The overwhelming 
majority of particles in the selected events are hadrons. These are known to suffer from hadronic interac- 
tions with the material in the detector. Thus a good description of the material in the detector is needed to 
get a good description of the track reconstruction efficiency. To quantify the influence of an imperfect de- 
scription of the detector description, in particular the material in the simulation, two different data-driven 
methods are used. The first reconstructs the invariant mass of mesons decaying to two charged pions; 
the second compares the track lengths in data and simulation. The /T*' mass method studies the mass 
as a function of the decay radius of the meson; it has greatest sensitivity to small radii, while the track 
length study probes the material description in the simulation in terms of nuclear interaction length {A) 
in the SCT detector. The combination of both methods provides good sensitivity throughout the silicon 
detectors. They allow us to constrain the material to better than 10% in the central barrel region and 
better than 30% at the highest I77I measured. The material uncertainty is the largest uncertainty in almost 
all regions of all distributions plotted in this paper. In the barrel region, the total uncertainty due to the 
material is 8% at low pT^, going down to 2% above 500 MeV. The uncertainty increases with increasing 
\r]\; the largest uncertainties are in the region 2.3 < I77I < 2.5: 15% in the first p-^ bin decreasing to 7% 
above 500 MeV. 

The track-fit probability cut has been found to offer powerful discrimination against tracks with 
mis-measured momenta. These are mostly very low momentum particles that are reconstructed with 
much higher momentum due to mis-alignment or nuclear interactions |3l Mis-measured tracks are seen 
predominantly at the edges of the 7/ acceptance where the distance between consecutive measurement 
points of the outer layer of the Pixel and the first layer of the SCT can reach up to ~ 1 m. The fraction of 
mis-measured tracks is observed to be significantly more in data than in Monte Carlo even after this cut is 
applied. Two different methods aie used to estimate the fraction of mis-measured tracks in data. The first 
compares the momentum obtained from the tracks reconstructed using only the SCT hit information with 
that obtained for fully reconstructed tracks. After normalising the number of well-measured tracks in MC 
to data, the scaling of the MC high-px tails needed to model the data is obtained. The second method 
uses the difference between data and MC seen in the tails of the do distributions at high pt; because mis- 
measured tracks tend to have poorly reconstructed d^. Again a scaling factor is obtained to scale the MC 
tails in order to describe the data. These two methods give very similar results. Both methods are used 
to obtain the systematic uncertainty for all but the outer-most regions in 77 where the effect is the most 
significant. In this region an additional method is used that compares the 77 distributions, normalised 
in the central region, in bins of /jt- The variation with pt; of the 77 distribution due to physics is small 
compared to the differences observed due to mis-measured tracks. The additional tracks at high I77I, high 
Pt: are considered to be due to mis-measured tracks and the fraction of mis-measured tracks in data is 
obtained. This third method gives the systematic uncertainty for the outer-most 77 bins. Averaged over 
the whole ?/ region, the fraction of mis-measured tracks in data is found to be negligible for px < 10 GeV, 
3% for 10 < /jt < 15 GeV and increases to 30% for 30 < /jt < 50 GeV. An additional systematic on 
the track reconstruction efficiency of 10% is taken for all tracks with pT^ > 10 GeV due to different 
efficiencies of the probability cut in data and MC. All systematic uncertainties on the mis-measured 
high-/7T tracks are taken as single-sided errors. 

Studies using Z ^ p/u events show that the resolution in data is about 10% worse than the nominal 
MC resolution above 10 GeV. The impact of a 10% Gaussian smearing of the reconstructed track pj in 

' Note that the momentum spectrum falls by many orders of magnitude in the measured range. 
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Systematic Uncertainty 


Size 


Region 


Material 


±2 - 15% 


decreases with px, increases with \t]\ 


prob. cut 


±10% 


flat, only for pj > 10 GeV 




±5% 


100 <pt < 150 MeV 


Resolution 


negligible 


0.15 <pj < 10 GeV 




-7% 


Pt > 10 GeV 


Track Selection 


±1% 


flat in pt and ij 


Truth Matching 


±1% 


only for = 2.36 TeV Pixel Tracks 


Efficiency correction factor 


±4% 


only for ^/s = 2.36 TeV ID Track 


Alignment and other high pj 


-3% to -30% 


only for pj > 10 GeV 
averaged over 77 , increases with increasing pj 



Table 5: The systematic uncertainties on the track reconstruction efficiency for = 0.9 TeV, - 
7 TeV and ^/s = 2.36 TeV Pixel Track and ID Track methods. Unless otherwise stated, the systematic 
is similar for all energies and phase-space regions. All uncertainties are quoted relative to the track 
reconstruction efficiency. 



MC is performed and found to have a 7% effect for the binning used in this paper. This effect is taken as 
a systematic uncertainty on tracks above 10 GeV. This systematic uncertainty is single-sided and added 
linearly with the systematic uncertainty due to the mis-measured high-px tracks. The effect on tracks 
below 10 GeV is found to be negligible. 

The Pj cut applied at various stages of the pattern recognition inside the track reconstruction algo- 
rithm introduces an inefficiency due to the momentum resolution. A different momentum resolution or a 
bias in the momentum estimation in data compared to MC can result in a change in the migration out of 
the first bin in pj (100 < pj < 150 MeV) and thus a gain or loss of observed tracks. The default migra- 
tion correction is derived using the resolution in Monte Carlo. The track pj resolution at the seed finding 
stage in Monte Carlo is increased by a very conservative 10 MeV, making the pj resolution effectively 
15 MeV instead of 10 MeV. The effect of this shift on the track reconstruction efficiency in the first pj 
bin is found to be about 5%; this difference is assigned as a systematic uncertainty. 

A detailed comparison of track properties in data and simulation is performed by varying the track 
selection criteria. The largest deviations between data and MC are observed at high rj and are found to 
be ~ 1%. For simplicity, a constant 1% uncertainty is assigned over ffie whole range. 

A summary of the track reconstruction systematic uncertainties is shown in Table [51 The total un- 
certainty due to the track reconstruction efficiency determination is obtained by adding all effects in 
quadrature except for tracks above 10 GeV where the resolution and mis-measured track effects are 
added linearly; asymmetric errors are considered for these effects. 

6.4 Track- Reconstruction Efficiency for the 2.36 TeV Data Sample 

Both the Pixel track and the ID track methods apply a data-driven correction to the primary track recon- 
struction efficiency, emc 

£{x) = SMcix) ■ Sconirj), (2) 

where emc is derived from nominal simulation at = 2.36 TeV. Here x is either boffi pj and ij for the 
ID track or only rj for the Pixel track method, as those are the parameters that the correction factors were 
found to depend on. 

The correction, ficorr. is derived from the reference dataset taken at ^/s = 0.9 TeV where the high 
voltage on the SCT was lowered for part of the run. 
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For the Pixel track method, gcorr is the ratio of the relative Pixel track reconstruction efficiency, Creh 
in data to simulation. The relative Pixel track efficiency is the efficiency to reconstruct a Pixel track if a 
track has been reconstructed using hits in the SCT and TRT detectors only. 

Figure shows the reiative Pixei track efficiency in data and simulation. The ratio of the two 
distributions, shown in the insert, is used to correct the track reconstruction efficiency for the Pixel track 
method at ^/s = 2.36 TeV. 

For the ID track method the efficiency derived from simulation with nominai conditions is corrected 
by ficorr to account for the iower SCT efficiency in standby mode. Figure |4j5 shows the distribution of 
the number of reconstructed tracks in data in both SCT configurations at = 0.9 TeV normalised to 
the same number of events satisfying the trigger requirement. The ratio of the number of reconstructed 
tracks with the SCT in standby, A^^'', to the number of reconstructed tracks with the SCT at nominal, 
N^°^, shown in the inset, is used to correct the track reconstruction efficiency for the ID track method at 

- 2.36 TeV: 

^coAri) - J^^^y (4) 
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Figure 4: Reiative efficiency of Pixel tracks in data and non-diffractive (ND) MC simulation at = 
0.9 TeV (a). Both Pixel track distributions are re-weighted to have the same beam spot distribution as 
the ^/s = 2.36 TeV data. The number of reconstructed ID tracks in data at = 0.9 TeV as a function 
of T] with the SCT in nominal and standby (b). The ID track distributions are normaiised to the number 
of events passing the trigger requirement. 



Systematic Uncertainties Most systematic uncertainties on the ID track reconstruction efficiency are 
similar to the full tracking at other energies. The major additional systematic uncertainty is due to the 
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efficiency correction factor for ttie SCT configuration. Tlie uncertainty due to tlie statisticai iimitations 
of ttie reference dataset is 2%. An additionai 3% uncertainty accounts for tlie extrapolation from = 
0.9 TeV to yfs = 2.36 TeV, which was estimated by compaiing the distributions of the number of ID 
tracks between = 0.9 TeV and = 2.36 TeV. The total uncertainty on the efficiency correction 
factor adds those two effects in quadrature to obtain a total uncertainty of 4%. 

The material uncertainty is estimated using a similar method as for the other energies; the absolute 
uncertainty is found to be 2% (3%) for the Pixel (ID) track reconstruction efficiency. The uncertainty 
is larger for ID tracks, because such tracks are sensitive to the material throughout the whole silicon 
detector. The uncertainty due to the momentum resolution is negligible because the phase-space cuts are 
sufficiently far from the track algorithm cuts. 

There is an additional 1% uncertainty on the Pixel track method due to the matching procedure. The 
relative Pixel track reconstruction efficiency differs from the primary efficiency due to material effects and 
contributions from non-primary tracks. There is an additional discrepancy of 4% in for 2.4 < \r]\ < 2.5 
that is assigned as a systematic uncertainty for those bins. At central rj the total uncertainty on the Pixel 
(ID) track reconstruction efficiency is estimated to be 3.4% (6%). Table[5]shows the track reconstruction 
systematics at ^/s = 2.36 TeV and the differences with respect to the uncertainties at other centre-of-mass 
energies are indicated. 



7 Correction Procedure 

The effect of events lost due to the trigger and vertex requirements is corrected using an event-by-event 
weight: 



"^ev^e?) 



1 1 



BSx _ 

£trig(«fel)'^vtx(«fel'^)' 



WtrkiPT, V) = • (1 - /nonp(PT)) " (1 - fokiiPJ, ?/)). 



where x is either the Az between tracks or the rj of the tracks, as described in Sec. 

The pj and rj distributions of selected tracks are corrected for using a track-by-track weight: 

1 
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where /nonp is the fraction of non-primary tracks determined as described in Sec. |5] 

The fraction of selected tracks passing the kinematic selection for which the corresponding primary 
particle is outside the kinematic range, fok^ipi, ??), originates from resolution effects and has been esti- 
mated from MC. The uncertainty on f^^ is mostly due to the resolution difference between data and MC. 
This uncertainty is negligible for all cases except at = 2.36 TeV for the Pixel track method where 
the uncertainty is estimated to be 1%, due to the poor momentum resolution of the Pixel tracks. No 
additional corrections are needed for the 77 distribution; the additional corrections needed for the other 
distributions are described in the following sections. 

For all distributions in all phase-space regions considered, closure tests are carried out. These are 
tests carried out on MC where the reconstructed samples are corrected according to the same procedure 
as used on the data; the resulting difference between the corrected distribution and the known particle 
level distribution is defined as the amount of non-closure; if the correction procedure were perfect, the 
non-closure would be zero. For this analysis, closure tests are carried out on all distributions in all 
phases-space regions and unless explicitly mentioned in the text the level of non-closure is less than 1%. 

7.1 Correction to ^ 

First, the observed Wsei distribution is corrected for the trigger and vertex reconstruction efficiencies. 
Then, an event-level correction is applied using Bayesian unfolding H6l to correct the observed track 
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multiplicity to the distribution of the number of primary charged particles, as follows. An unfolding 
matrix, Mch,sei» is defined that expresses the probability that a given selected track multiplicity, after all 
other event-level corrections are applied, Wseh is due to n^h primary particles. This matrix is normahsed 
such that the number of events does not change except for the rare cases where n^ei > rich and n^h is 
below our acceptance selection. This matrix is populated from MC09 MC and appUed to data to obtain 
the observed «ch distribution. The resulting distribution is then used to re-populate the matrix and the 
correction is re-appUed. This procedure is repeated without a regularisation term and converges after 
four iterations in data; convergence is defined as the first iteration in which thex^ difference between the 
result of the unfolding and the input distribution for that iteration is less than the number of bins used in 
the unfolding. 

After the n^ei distribution has been unfolded, the resulting charged particle multiplicity distribution 
is corrected for events migrating out of the selected kinematic range (rich > X), which the matrix does 
not account for. This is achieved by adding an additional term to the correction. The correction terms for 
the phase-space regions with nch > 2 is 

1 /(I - (1 - etrk)"* - nch • etrk ■ (1 - s^yd^"^^-'^) (5) 

where etrk is the mean effective track reconstruction efficiency for a given Uch bin. Corresponding terms 
are used for the other phase-space regions. This track reconstruction efficiency can in principle be differ- 
ent for each Uch bin, but the difference is found to be small and thus the mean effective track reconstruction 
efficiency for lowest rich bin is used. 

Systematic Uncertainties The systematic uncertainties on the unfolding procedure are obtained by 
modifying the input distributions as described below, applying the unfolding procedure and comparing 
the output to that obtained when using the nominal input; the matrix and the correction factors are not 
modified. 

There are two sources of systematic uncertainties considered. One of them is due to the track re- 
construction efficiency uncertainties while the second one accounts for the different spectrum recon- 
structed in data and MC. The first source of uncertainty is estimated by starting from the observed Wgei 
spectrum in data; tracks are randomly removed from the distribution according to the mean pj and rj of 
the tracks for each value of nsei and the uncertainty on the track reconstruction efficiency for those pj and 
77 values. A new input distribution is obtained, put through the unfolding procedure and the difference 
with respect to the nominal nch distribution is taken as a systematic uncertainty. The uncertainty is then 
symmetrised. The uncertainty on «ch due to the uncertainty on the track reconstruction efficiency is found 
to be ~ 3% to ~ 25% at ^/s -1 TeV in the most inclusive phase-space region, nch > 1, pi > 100 MeV, 
|77| < 2.5. 

The other source of uncertainty originates from the unfolding method that is carried out in a single 
dimension at a time, in this case «ch- There is some dependency on the spectrum of the MC sample 
used to populate the matrix, due to the strong dependence of the track reconstruction efficiency on pj. To 
investigate this effect, the average track reconstruction efficiency derived using the pi spectrum in data 
and that obtained from MC are compared. The difference in these two mean efficiencies is then treated in 
the same way as the uncertainty on track reconstruction efficiency, described in the previous paragraph. 
This uncertainty is taken as being asymmetric; only the contribution from a shift of the spectrum in the 
direction of the data is taken. The mean value is kept as that given by the nominal px spectrum in MC. 
The uncertainty varies with increasing nch from -2% to -1-40% at yfs - 1 TeV in the most inclusive 
phase-space region. 

The only additional systematic uncertainty due to the tuning of the track reconstruction efficiency 
is due to the difference between the bias introduced by the vertex correction in MC and data. The 
estimation of this error is done by comparing the Azo distribution in n??=2 between data and MC. The 



21 



Azo distribution is a very good probe of the correlation between ?isei/'^ch and n^^^ as events with high ?isei 
tend to have small Azo values while events with ?isei<2 tend to have large Azo- Very good agreement 
is found between data and MC. Re- weighting the Azo distribution in MC to match data or applying 
the vertex correction extracted from data to the MC closure test leads to a systematic uncertainty of 
the order of 0.1% for ?ich=2 where this effect is most pronounced. As this error is much smaller than 
other systematic uncertainties considered, it is neglected. The systematic uncertainty due to track-track 
correlation in a single event is small and is neglected everywhere in this analysis. 

7.2 Corrections to A^ev 

The total number of events, A^ev> used to normalise the final distributions, is defined as the integral of the 
rich distributions, after all corrections are appUed. 

Systematic Uncertainties The systematic uncertainties on A^ev are obtained in the same way as for the 
«ch distributions. Only those systematics affecting the events entering or leaving the phase-space region 
have an impact on A^ev The total uncertainty on A^ev at = 7 TeV for the most inclusive phase-space 
region is 0.3%, due mostly to the track reconstruction efficiency. At = 2.36 TeV the total uncertainty 
on A^ev is 1.4% for the Pixel track and 2.6% for the ID track methods. 

7.3 Corrections to - ^ 

The tracks are first corrected for the event level inefficiencies of the trigger and the vertex reconstruction. 
Then the tracks are corrected for the track reconstruction inefficiencies, non-primary track contamination 
and out of kinematic range factors. Finally, a similar unfolding method to that used on the rich distribution 
is used to correct the measured track pj to the primary particle momentum. More bins are used for the 
unfolding than are shown in the final distributions; this is necessary in order to avoid amplification of 
small data MC differences with successive iterations, causing large fluctuations. For this distribution four 
iterations are required before convergence is reached; convergence is defined as for the rich distribution. 

Systematic Uncertainties In order to estimate the effect on the final pj distributions of the uncertain- 
ties affecting the correction steps prior to the unfolding, the unfolding procedure is re-run on the corrected 
Pj distribution shifting the distribution used as input to the unfolding procedure by the systematic uncer- 
tainties. This new pj distribution is put through the unfolding procedure and the difference with respect 
to the nominal corrected pj spectrum is taken as a systematic uncertainty. 

The high-/7T systematic uncertainties are obtained using the MC samples. The systematic uncertainty 
associated to the mis-measured high-px tracks is obtained by scaling the number of mis-measured tracks 
in MC to match those found in data. This new input distribution is put through the unfolding procedure 
and the final difference with respect to the nominal MC is taken as a systematic uncertainty. The system- 
atic uncertainty associated to the resolution is obtained by smearing the well-measured tracks, in MC, by 
the resolution uncertainty obtained in Sec. 16.31 The effect on the final unfolded distribution is taken as a 
systematic uncertainty. Those two high-jox systematics are added linearly. Both cause only single-sided 
variations. This combined uncertainty is measured to be from -10% for pj = 10 GeV to -30% for the 
last Pj bin (30 < pj < 50 GeV) at = 7 TeV for the rich > 2, pj > 100 MeV phase-space region. 
The variations for other phase-space regions at this energy are similar. At = 0.9 TeV this uncertainty 
is found to be -20% for all three bins above /jt of 10 GeV. 

In order to assess the stability of the results under varying starting hypotheses for the MC spectrum 
used to fill the matrix, a flat initial prior is used as an input. While convergence is only typically reached 
after seven iterations, instead of three for the nominal prior, the final difference in the unfolded spectra is 
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small. The difference between the resulting distribution obtained with a flat prior and that obtained with 
the MC pj spectrum as a prior is taken as a systematic uncertainty. At ^/s = 7 TeV this uncertainty is 
less than 2% for nearly all pT bins, with the exception of a couple of bins around changes in bin width, 
where the effect is 3-5%. At = 0.9 TeV, due to more hmited statistics in the MC, the largest change 
seen is 7% with a few others around 3-4%. 

7.4 Mean pj versus rich 

The correction procedure for the (pj) vs. rich distribution is designed to correct separately two compo- 
nents: Yji Pjii) vs. rich and X,- 1 vs. rich and take the ratio only after all corrections are applied. The sum 
is over all tracks and all events; the first sum is the total pj of all tracks in that bin in rich', the second 
sum represents the total number of tracks in that bin. The sums will be referred to as the numerator and 
denominator, respectively. Each of these distributions, li/PxCO and 1> is corrected in two steps. 

First the two distributions as a function of ?isei are corrected on a track-by-track basis by applying the 
appropriate track weights; this track-by-track correction is applied to the data distribution and thus no 
longer relies on the pj spectrum of the MC. Second, the matrix obtained after the final iteration of the 
rich unfolding described in Sec. 17. H is applied to each of the distributions to unfold «sei to rich- Finally, the 
ratio of the two distributions is taken to obtain the corrected (pj) vs. rich distribution. For this distribution 
we exclude tracks with pj > V^/2 as they are clearly un-physical; this removes 1 track at = 0.9 TeV 
and 1 track at = 7 TeV. 

This unfolding procedure assumes that the tracking efficiency depends only on pj and rj and is 
independent of the track particle multiplicity, and that the pj spectrum of the tracks in events that migrate 
back from a given ?isei bin to a given «ch bin is the same as the pj spectrum of tracks in events in the 
corresponding Wsei bin. The fact that these assumptions are not completely valid is taken as a systematic 
uncertainty. This uncertainty is obtained by looking at the non-closure of the corrected distribution in 
the MC. This residual non-closure is, we believe, a consequence of the two main assumptions. A full 
parametrisation of the track reconstruction efficiency in terms of pj, rj and «ch would remove the need for 
the first assumption, while a full two-dimensional unfolding as a single step where the two dimensions 
were pj and rich would remove the need for the second. Both of these are beyond the scope of the current 
paper. In order to understand if the amount of non-closure is a realistic estimate of the uncertainty on 
the method when applied to data, in particular to investigate its dependence on the pj spectrum, the 
whole unfolding procedure is carried out using pythia6 DW tune samples and the pythiaS samples; we 
varied both the input distribution and the matrix used to do the unfolding. The level of non-closure is 
found to be similar to that obtained with the MC09 pythia6 samples. We thus conclude that the level of 
non-closure is not strongly dependent on the pj spectrum. This allows us to use the residual non-closure 
as a systematic uncertainty on the unfolding method as described in the next section. 

Systematic and Statistical Uncertainties For the calculation of the statistical uncertainty, the full 
correlation between the tracks inside the same event was not computed. The statistical uncertainty in the 
numerator and denominator are computed separately then added in quadrature after taking the ratio. This 
is found to be a conservative estimate of the uncertainty. 

Systematic uncertainties considered for the (pt) vs. rich distribution are either due to assumptions 
made during the correction procedure or to uncertainties on quantities taken from the MC and used during 
the correction procedure. 

The first category refers to the assumptions on the method, the effects of which are visible in the 
closure test. To account for these imperfections, we apply a systematic uncertainty of 2%, which covers 
the non-closure in MC, except for the highest rich bin and the first few rich bins in some of the phase- 
space regions. For these cases a larger systematic uncertainty is applied to cover the non-closure. For the 
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analyses with pj > 500 MeV, where the size of a non-closure is larger, a 3% systematic error is applied 
in the ?ich=l bin. This systematic uncertainty also covers the difference in the non-closure between 
samples created using MC09 (default) and those with DW tune of pythia6 and pythiaS. In the correction 
procedure we use the approximation that «sei= «fgp The effect of such an approximation is studied on 
simulation and found to be negligible with respect to the other sources of uncertainty. 

The second category comprises uncertainties on the track correction weights Ws^,{f^^^^) and w^^ipi, rf) 
and on the migration probabilities obtained from the unfolding matrix. The dominant systematic uncer- 
tainties that affects both the track corrections weights and the migration probabilities are the same as 
those affecting the n^h distribution unfolding: the uncertainty on the track reconstruction efficiency and 
the effect of the difference in the /jj spectra between data and MC. These uncertainties are propagated 
by varying the input distribution for both the Yji Piii) vs. «sei and X,- 1 vs. ?isei. 

Smaller effects are also studied, for example the uncertainty on the rate of non-primary tracks and 
the effect of the systematic uncertainties affecting the high-pj tracks mentioned in Sec. 16.31 Excluding 
the systematic uncertainties due to the assumptions made during the correction procedure, the systematic 
uncertainties are between 0.5% and 2% for all bins in rich, all energies and all phase-space regions. 

7.5 Correction for Different Minimum n^h Requirements 

The only difference in the correction procedure from track to particle level for «ch ^ 6 with respect to 
rich > 1 is the need for an additional correction that takes into account the effect on the tracks due to the 
tighter cut on both the number of tracks and number of particles. 

The rich distribution and the number of events A^ev are obtained by correcting and unfolding the mul- 
tiplicity distribution of the whole spectrum and then applying the higher rich cut on the final distribution. 
For the px and rj track distributions an extra correction is needed. For events with Wsei ^ 6, the tracks 
are added to the distribution as for all other phase-space regions; a weight corresponding to the product 
of the track (wtrk) and event weights (ifeu) is applied. For events with ?isei < 6 the tracks are added to the 
distribution with an additional weighting factor, u;„^^<6 that represents the probability that a track from 
an event with ^sei tracks is from an event with rich ^ 6. This additional weight is taken from the final 
rich unfolding matrix, after the final iteration; each column in the matrix represents the probability that 
an event with rigci tracks has rich particles. The total probability {p{nch > 6 | «sei)) for a given ?isei < 6 
is therfor the sum over the matrix elements for rich ^ 6 

where Mn^^ n^^j is the entry in the unfolding matrix for rich and n^ci- This weight is about 65% for ?isei = 5 
and rapidly drops to 1% for «sei - 2. 

Systematic Uncertainties All uncertainties related to the distributions with the lower rich cut are taken 
into account. In addition, an extra systematic uncertainty due to the uncertainty on the track reconstruc- 
tion efficiency is needed for the correction to higher rich selection. By varying the track reconstruction 
efficiency down by its uncertainty, different w,i^^<6 weights are obtained. The shift in the resulting rich 
distribution is symmetrised and taken as an additional systematic uncertainty. 

7.6 Extrapolation to pr = 

Comparing the results in our well-defined phase-space regions to other inclusive measurements from 
other experiments requires additional model-dependent corrections. One such correction is described 
here, but applied only for comparative purposes. This particular correction is derived to extrapolate 
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the average multiplicity in the phase-space region with the lowest measured pj to the multiplicity for 
all pt > 0. No attempt is made to correct for the rich ^ 2 requirement. Results are quoted for the 
average multiplicity in the rapidity interval \t]\ < 2.5 and are not considered to be the main results of this 
paper. This correction is obtained using three independent methods: fitting the pj spectrum to a given 
functional form, assuming a flat distribution at low pj in the observed fully corrected ^ • distribution 
and obtaining the correction factor from the AMBTl pythia6 MC. 

In the first method, the corrected pj spectrum is fit with a two-component Tsallis II47II48II distribution 
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where mripT) is the transverse mass mr - yjpj + and mo is the particle rest mass mo - {m„,mp} 
and dNcb/dy\y=oj, Tt and are the six parameters of the fit. ?/' represents the pseudorapidity at the edge 
of our acceptance, 77 - 2.5. dNch/dy\y=o represents the integrated yield of the particle production at 
mid-rapidity, but is left here as a free parameter of the fit. Mesons (pions and kaons) are merged into 
a single Tsallis function since there is insufficient information in the measured distribution to fit three 
independent shapes. The tanh"' factor accounts for the variation in E/p of each track over the entire 
measured pseudorapidity range. It is derived by integrating ^drj over \r]\ < 2.5. 

From this functional form and using the parameters obtained from the fit, the fraction of particles 
with Pi < 100 MeV is extracted. This procedure gives the correction factor to be applied to the mean 
charged-particle multiplicity per unit ij , averaged over \r]\ < 2.5, in order to get the inclusive multiplicity. 
The correction factor from pj > 100 MeV to pj > MeV is found to be 1.065 at = 0.9 TeV and 
1.063 at - 7 TeV. 

The second method assumes that the — • distribution is flat at low pj. One can thus use the value 
of this distribution in the lowest pj bin (100 < pj < 150 MeV) to extract the value for tracks below 
100 MeV. From this assumption, the fraction of particles below 100 MeV and the scale factor used to 
correct our observed distributions are derived. The scale factors are found to be 1 .068 at = 0.9 TeV 
and 1 .065 at = 7 TeV. The third and final method simply obtains the correction factor using one of the 
MC models. AMBTl PYTmA6 is chosen; the correction factors are found to be 1.055 at = 0.9 TeV 
and 1.051 at = 7 TeV. We chose to use the scale factor obtained from the functional form fit as 
the central value and consider the difference between this and the other two methods as a systematic 
uncertainty. 



Systematic Uncertainties Several sources of systematic uncertainty on the calculated scale factor are 
considered. The dominant uncertainty comes from the difference in the scale factors obtained from the 
three different extrapolation methods. The largest difference between the value obtained from the fit 
and the values from the MC and from the flat extrapolation is considered as the uncertainty and then 
symmetrised. This uncertainty is found to be 0.007 at = 0.9 TeV and 0.012 at ^/s = 7 TeV. 

The other sources of uncertainty are related to the fitting procedure such as the variation within the 
uncertainty on the fit parameters and the variation due to a change of the the fit range. All sources of 
uncertainty are assumed to be uncorrelated and thus added in quadrature. The final scale factors, with 
total uncertainty, are then 1.063 ± O.OIV at = 7 TeV and 1.065 ± 0.01 Itot at = 0.9 TeV. 
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8 Total Systematic Uncertainties 



The individual sources of systematic uncertainties have already been discussed in previous sections. 
The effect on the final distribution from each source is treated independently and propagated to the 
final distributions; the total error is the sum in quadrature from the different sources, unless explicitly 
mentioned in the text. In most bins of all distributions the largest uncertainty comes from the track 
reconstruction efficiency. The uncertainties at = 2.36 TeV are larger than at the other two energies 
due to the uncertainties related to the operation of the SCT at reduced bias voltage during 2.36 TeV data 
taking. The total uncertainties are shown as shaded bands in the final distributions presented in the next 
section. 

9 Results and Discussion 

The corrected distributions for primary charged particles for events in three separate phase-space regions 
are shown in Fig. [5] to \T3\ The results are compared to predictions of models tuned to a wide range 
of measurements. The measured distributions are presented as inclusive-inelastic distributions within a 
given phase-space region with minimal model-dependent corrections to facilitate the comparison with 
models. 

9.1 Charged-Particle Multiplicities as a Function of the Pseudorapidity 

Figures [5] and [6] show the charged-particle multiplicity as a function of pseudorapidity. Figure [5] shows 

the distribution at all three centre-of-mass energies in the phase-space region, ?ich ^ I.Pt > 500 Me V, 
I77I < 2.5. The mean particle density is roughly constant for \t]\ < 1.0 and decreases at higher values 
of It/I. There is little shape variation between the models except for the DW pythia6 tune which has a 
flatter spectrum and a more pronounced dip at central 77 , especially at low V^. At all three energies 
the AMBTl pythia6 tune gives the best shape and normalisation description of the data, although it was 
tuned for rich > 

Figure [6^ and b show the rj distributions for the most inclusive phase-space region, «ch ^ 2, 
Pt > 100 MeV, I77I < 2.5. There is less rj variation than in the previous figure. At 900 GeV there 
is very little difference between the models both in shape and normalisation with the exception of phojet 
which shows an excellent agreement with the data; the other models show on average too few particles. 
The shape of the distribution is reasonably well described by all models. At 7 TeV again the shapes seem 
to all model reasonably well the observed spectrum, but at this energy the difference in normalisation 
among the models varies more widely and no model reproduces the data. 

Figure |6]: and d show the rj distributions for the phase-space region with the least amount of diffrac- 
tion, rich ^ 6, ;?T > 500 MeV, \r]\ < 2.5. The distributions in this phase-space region have the largest 
drop at high |?7|. All but pythia6 DW and phojet at ^/s = 1 TeV show reasonable agreement in both shape 
and normalisation at both energies. 

9.2 Charged-Particle Multiplicities as a Function of the Transverse Momentum 

Figures |7] and [8] show the charged-particle multiplicities as a function of the transverse momentum. 
The first of these figures shows all three centre-of-mass energies considered in the phase-space region 
'^ch ^ ^, Pt > 500 MeV and |?7| < 2.5. The observed pj spectrum is not described by any of the 
models over the whole range. The region that the models have the most difficulty describing is the region 
above 1 GeV. 

Figures [8^ and b show the charged-particle multiplicities in the most-inclusive phase-space region. 
At 900 GeV phojet describes the data best over the whole range even though the agreement is still not 
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excellent. The other models tend to under-predict the number of low pi particles while at higher pj the 
models vary widely. At 7 TeV the effect at low pj is more pronounced, while at high pj the agreement 
of pythiaS and phojet with the data is quite good. The AMBTl and MC09 tunes of pythia6 predict too 
many particles at higher pj. 

Figures [St and d show the charged-particle multiplicities with the smallest contribution from diffrac- 
tive events. This distribution carried the most weight in the AMBTl tune. Considerable improvement in 
the agreement with data is seen between the older MC09 and the newly tuned AMBTl but the parameters 
varied in this tune were not sufficient to describe the full spectrum. 

9.3 Charged-Particle Multiplicity Distribution 

Figure|y]shows the charged-particle multiplicity distributions for ?ich ^ ^,Pt > 500 MeV and I?/! < 2.5 
at all three centre-of-mass energies. At low number of charged particles, all models predict more events 
than observed in data, which is compensated by an under-prediction in the tails of the distributions. 
It should be noted that due to the normalisation, 1/A^ev> a deviation observed in one region needs to 
be compensated for by one in the other direction somewhere else. Although the predictions of phojet at 
0.9 TeV model the data reasonably well, at 2.36 TeV and 7 TeV they do not model the observed spectrum. 
The new AMBTl pythia6 tune seems to provide the best agreement with data. 

Figures [TOb and b show the distribution for the most inclusive phase-space region. Here the variations 
between models at both low and high values of rich are increased and no model predicts the observed 
spectra. 

Figures [TOb and d show the distribution for the diffraction-reduced phase-space region. The distribu- 
tions are very similar to those in Fig. |9] with a cut at rich ^ 6; only the normalisation is different between 
the plots. The errors are also recomputed as there is a larger cancellation between the numerator and 
denominator for this phase-space region. 

9.4 Average Transverse Momentum as a Function of the Number of Charged Particles 

The final set of distributions discussed in the main part of this paper is the average transverse momentum 
as a function of particle multiplicity. The measurement of (pj) as a function of charged multiplicity at 
= 2.36 TeV is not shown because different track reconstruction methods are used for determining the 
Pj and multiplicity distributions, as discussed in Sec. 14.3.2] Figure [TT] shows the results for events with 
«ch ^ 1> Pt > 500 MeV and \t]\ < 2.5. At 900 GeV the slope vs. «ch for high values of «ch seems to 
be well described by most models but the absolute value is best modelled by pythia6 DW. At the highest 
centre-of-mass energy above 20 particles the models vary widely both in slope and in absolute value; 
at low values of «ch none of the models describe the data very well. In the more inclusive phase-space 
region, Fig.[T2k and b, the models vary widely, especially at high V^. 

9.5 drichldrj at 77 = 

The mean number of charged particles in the central region is computed by averaging over I;;! < 0.2. 
The values for all three phase-space regions and all energies available are shown in Fig.[T3]and in Tabled 
The result quoted ni ^fs = 2.36 TeV is the value obtained using the Pixel track method. The phase-space 
region with largest minimum pj and highest minimum multiplicity {pi > 500 MeV; rich ^ 6), which 
is the region with the least amount of diffraction, is the one where the models vary the least and the 
energy extrapolations of most models agree the best with the data. However, in this region the energy 
extrapolation of pythia6 and phojet do not agree with the data. For the most inclusive measurements, 
none of the models agree with the data and the spread at 7 TeV in the expected values is almost one third 
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of the mean predicted value. The observed value is significantly higher at this energy than any of the 
models. 



Phase-Space Region 


Energy 
(TeV) 


dnch/dr] at rj - 
Measured 




PYTHIA6 AMBTl MC 


rich > 2,pT > 100 MeV 


0.9 

7 


3.483 ± 0.009 (Stat) ± 0.106 (syst) 
5.630 ± 0.003 (Stat) ± 0.169 (syst) 


3.01 
4.93 


rich > ^,Pt > 500 MeV 


0.9 
2.36 
7 


1.343 ± 0.004 (Stat) ± 0.027 (syst) 
1.74 ± 0.019 (Stat) ± 0.058 (syst) 
2.423 ± 0.001 (Stat) ± 0.050 (syst) 


1.28 
1.70 
2.36 


rich > 6, pt > 500 MeV 


0.9 

7 


2.380 ± 0.009 (Stat) ± 0.027 (syst) 
3.647 ± 0.002 (Stat) ± 0.052 (syst) 


2.33 
3.63 



Table 6: drich/dr] at j] - for the three dilferent phase-space regions considered in this paper for the 
energies where results are available. For MC, sufficient statistics were generated such that the statistical 
uncertainty is smaller than the last digit quoted. 



9.6 Extrapolation to pt = 

The mean multiplicities of charged-particles with pj > 100 MeV within the full \r]\ < 2.5 region are 
computed as the mean of the distributions shown in Fig. [6^ and b. They are found to be 3.614 ± 0.006 
(Stat) ± 0.170 (syst) at ^ 0.9 TeV and 5.881 ± 0.002 (stat) ± 0.276 (syst) at = 7 TeV. Multiplying 
these numbers by the model-dependent scale factors obtained in Sec. l7.6[ the averaged inclusive charged- 
particle multiplicity for events with two or more particles is found to be 3.849 ± 0.006 (stat) ± 0. 185 (syst) 
at ^ 0.9 TeV and 6.252 ± 0.002 (stat) ± 0.304 (syst) at yfs ^ 7 TeV. This result is interpreted as 
the average total inelastic multiplicity for events with two or more particles within I77I < 2.5. Figure [14] 
compares these results to recently published ALICE results f5"61 for inclusive inelastic as well as inelastic 
with more than one particle. The ALICE results are quoted as averages over |7/| < 1.0and|?7| < 0.5, 
respectively. 

10 Conclusions 

Charged-particle multiplicity measurements with the ATLAS detector using the first collisions delivered 
by the LHC during 2009 and 2010 are presented. Based on over three hundred thousand proton-proton 
inelastic interactions at 900 GeV, just under six thousand at 2.36 TeV and over ten million at 7 TeV, the 
properties of events in three well-defined phase-space regions were studied. The data were corrected 
with minimal model dependence to obtain inclusive distributions. The selected kinematic range and 
the precision of this analysis highlight clear differences between Monte Carlo models and the measured 
distributions. In all the kinematic regions considered, the particle multiplicities are higher than predicted 
by the Monte Carlo models. 

The three different phase-space regions studied, from the most inclusive to the one with the smallest 
diffiactive contribution, highlight various aspects of the charged-particle spectra. In general, the agree- 
ment between the models and the data is better in the phase-space regions with higher minimum pj 
cutoff, where diffractive contributions are less significant. 

For the = 0.9 TeV measurements with the pj threshold of 500 MeV, these results supersede the 
results presented in HI. 
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Figure 5: Charged-particle multiplicities as a function of the pseudorapidity for events with n^h > 1, 
Pt > 500 MeV and < 2.5 at - 0.9 TeV (a), - 2.36 TeV (b) and = 7 TeV (c). The dots 
represent the data and the curves the predictions from different MC models. The vertical bars represent 
the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties added in 
quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio histograms 
refer to the bin centroids. 
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Figure 6: Charged-particle multiplicities as a function of the pseudorapidity for events with rich ^ 2, 
Pt > 100 MeV (a,b) and n^h > 6, pj > 500 MeV (c,d) and \r]\ < 2.5 at = 0.9 TeV (a,c) 
and ^fs - 1 TeV (b,d). The dots represent the data and the curves the predictions from different MC 
models. The vertical bars represent the statistical uncertainties, while the shaded areas show statistical 
and systematic uncertainties added in quadrature. The bottom inserts show the ratio of the MC over the 
data. The values of the ratio histograms refer to the bin centroids. 
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Figure 7: Charged-particle multiplicities as a function of the transverse momentum for events with 
«ch > ^, Pt > 500 MeV and < 2.5 at - 0.9 TeV(a), ^/s = 2.36 TeV(b) and - 7 TeV(c). 
The dots represent the data and the curves the predictions from different MC models. The vertical bars 
represent the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties 
added in quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio 
histograms refer to the bin centroids. 
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Figure 8: Charged-particle multiplicities as a function of the transverse momentum for events with 
rich > 2, pT > 100 MeV (a,b) and n^h > 6, pj > 500 MeV (c,d) and \t]\ < 2.5 at = 0.9 TeV 
(a,c) and ^[s -1 TeV (b,d). The dots represent the data and the curves the predictions from dilferent MC 
models. The vertical bars represent the statistical uncertainties, while the shaded areas show statistical 
and systematic uncertainties added in quadrature. The bottom inserts show the ratio of the MC over the 
data. The values of the ratio histograms refer to the bin centroids. 
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Figure 9: Charged-particle multiplicity distributions for events with rich ^ ^, Pt > 500 MeV and 
\r]\ < 2.5 at ^/s ^ 0.9 TeV(a), = 2.36 TeV(b) and = 7 TeV(c). The dots represent the data 
and the curves the predictions from different MC models. The vertical bars represent the statistical 
uncertainties, while the shaded areas show statistical and systematic uncertainties added in quadrature. 
The bottom inserts show the ratio of the MC over the data. The values of the ratio histograms refer to the 
bin centroids. 
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Figure 10: Charged-particle multiplicity distributions for events with rich ^ ^, Pt > 100 MeV (a,b) 
andnch > 6, pj > 500 MeV (c,d) and |?7| < 2.5 at ^ 0.9 TeV (a,c) and ^ 7 TeV (b,d). The dots 
represent the data and the curves the predictions from different MC models. The vertical bars represent 
the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties added in 
quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio histograms 
refer to the bin centroids. 
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Figure 1 1 : Average transverse momentum as a function of the number of charged particles in the event 
for events with Wch > 1, Pt > 500 MeV and < 2.5 at ^ 0.9 TeV(a), and ^ 7 TeV(b). 
The dots represent the data and the curves the predictions from different MC models. The vertical bars 
represent the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties 
added in quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio 
histograms refer to the bin centroids. 
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Figure 12: Average transverse momentum as a function of the number of charged particles in the event 
for events with «ch > 2, pT > 100 MeV and \r]\ < 2.5 at a/s ^ 0.9 TeV (a) and ^ 7 TeV (b). 
The dots represent the data and the curves the predictions from different MC models. The vertical bars 
represent the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties 
added in quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio 
histograms refer to the bin centroids. 
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Figure 13: The average charged-particle multiplicity per unit of rapidity for 77 = as a function of the 
centre-of-mass energy. The results with rich ^ 2 within the kinematic range pj > 100 MeV and 
\t]\ < 2.5 are shown alongside the results with «ch ^ 1 within the kinematic range pj > 500 MeV and 
I77I < 2.5 at 0.9, 2.36 and 7 TeV. The data are compared to various particle level MC predictions. The 
vertical error bars on the data represent the total uncertainty. 
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Figure 14: The average charged-particle multiplicity per unit of rapidity as a function of the centre-of- 
mass energy. The ATLAS results are for «ch ^ 2 in the region I77I < 2.5. For comparison ALICE results 
for «ch ^ 1 in the region |?7| < 1.0and?ich ^ in the region I77I < 0.5 are shown. It should be noted 
that the ALICE points have been slightly shifted horizontally for clarity. The data points are compared 
to pythia6 AMBTl predictions for the same phase-space regions. 
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A Distributions Used in AMBTl Timing 



Table |7] and [8] show the list of all distributions from ATLAS and the Tevatron, respectively, used in tl 
ATLAS Minimum Bias Tune 1 (AMBTl). The Analysis column refers to the event selection used 
the particular analysis. The Tuning range column refers to the portion of the phase-space region that 
considered for the tune. 



Analysis 



Observable 
-m-;r~ 



Tuning range 



1 

J_ . _}_ . d-jVch 

J_ . ^ 



ATLAS 0.9 TeV, minimum bias, «ch ^ 6 
ATLAS 0.9 TeV, minimum bias, «ch ^ 6 
ATLAS 0.9 TeV, minimum bias, rich > 6 
ATLAS 0.9 TeV, minimum bias, rich ^ 6 



drjdpj 
Ich 

<Pt) vs. rich 

(f^) vs. p'^'^'' (towards) 
< W> ^t"' (transverse) 

<Tif > ^t'" (towards) 
^4^^ ^t""^ (transverse) 

i M7h 



-2.5 < ?7 < 2.5 
Pt>5.0 GeV 
rich > 20 
rich > 10 



ATLAS 0.9 TeV, UE in minimum bias 
ATLAS 0.9 TeV, UE in minimum bias 
ATLAS 0.9 TeV, UE in minimum bias 
ATLAS 0.9 TeV, UE in minimum bias 
ATLAS 0.9 TeV, UE in minimum bias 
ATLAS 0.9 TeV, UE in minimum bias 



p}^^'^ > 5.5 GeV 
/jJp^'^ > 5.5 GeV 
pl^'"^ > 5.5 GeV 
pJ^^^'^ > 5.5 GeV 



A'cv d;; 
J_ . -L- . d^jVch 
Ncv Inpf drjdpj 

1 . dWev 

N„ dn^h 

(Pt) vs. rich 



Pi 
p}^ 



> 5.5 GeV 



> 5.5 GeV 



ATLAS 7 TeV, minimum bias, rich ^ 6 
ATLAS 7 TeV, minimum bias, rich ^ 6 
ATLAS 7 TeV, minimum bias, rich ^ 6 
ATLAS 7 TeV, minimum bias, «ch ^ 6 



1 / "cn 

< W> ^t'" (towards) 
< W> ^t"* (transverse) 

<w>^^- /^T^^away) 

<W > ^t'" (towards) 
(^i^^ '^^^ ^t""^ (transverse) 



-2.5 < ?7 < 2.5 
Pt>5.0 GeV 
Wch > 40 
rich > 10 



ATLAS 7 TeV, UE in minimum bias 
ATLAS 7 r^y, UE in minimum bias 
ATLAS 7 TeV, UE in minimum bias 
ATLAS 7 TeV, UE in minimum bias 
ATLAS 7 TeV, UE in minimum bias 
ATLAS 7 TeV, UE in minimum bias 



di^d<f) ' ' 1 ^ 

vs. p'^-^ (away) 



„lead 
Pj 


> 


10 GeV 


„lead 
Pj 


> 


10 GeV 


p'r" 


> 


\OGeV 


„lead 
Pj 


> 


10 GeV 


p'r' 


> 


10 GeV 


„lead 
Pj 


> 


10 GeV 



Table 7: ATLAS observables and ranges of distributions used in the AMBTl tuning. 
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Observables 

CDF Run I underlying event in dijet events H9l (leading jet analysis) 

A^ch density vs. leading jet pj (transverse), JET20 

A^ch density vs. leading jet pj (toward), JET20 

A^ch density vs. leading jet pj (away), JET20 

Yj Pt density vs. leading jet pj (transverse), JET20 

X Pt density vs. leading jet pj (toward), JET20 

X Pj density vs. leading jet pj (away), JET20 

A^ch density vs. leading jet pj (transverse), min bias 

A^ch density vs. leading jet pj (toward), min bias 

A^ch density vs. leading jet pj (away), min bias 

Pt density vs. leading jet pj (transverse), min bias 

Pt density vs. leading jet pT (toward), min bias 
X Pt density vs. leading jet pT (away), min bias 
Pt distribution (transverse), leading pt > 5 GeV 
Pt distribution (transverse), leading pT > 30 GeV 

CDF Run I underlying event in MIN/MAX-cones OHl ("MIN-MAX" analysis) 

(/jj?^"} vs. = 1800 GeV 

(p^^) vs. E^^^'^, - 1800 GeV 

(pf^) vs. E^^^'^, = 1800 GeV 

<A^max> vs. £■5,^^'', - 1800 GeV 

<A^min> vs. E^^^'^, ^f] = 1800 GeV 

Swiss Cheese p'^"" vs. E^^""^ (2 jets), = 1800 GeV 

(;7™> vs. = 630 GeV 

<p™"> vs. E^^'"^, - 630 GeV 

(pf^) vs. E^^""^, VI = 630 GeV 

Swiss Cheese p'^'^ vs. g'^"'^ (2 jets), = 630 GeV 

DO Run II dijet angular correlations 115111 
Dijet azimuthal angle, /j^'"' e [75, 100] GeV 
Dijet azimuthal angle, 7?^^" € [100, 130] GeV 
Dijet azimuthal angle, 7?^^" € [130, 180] GeV 

Dijet azimuthal angle, /p^^" >180 GeV 

CDF Run II minimum bias [|52l 
(Pt) of chai-ged particles vs. A^ch> = I960 GeV 
CDFRunlZpT [53] 
, = 1800 GeV 



Table 8: Tevatron datasets used in the AMBTl tuning. No specific cuts on the tuning ranges were made. 
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B Additional Phase-Space Regions 



Two additional phase-space regions are considered in this appendix: 

• at least twenty charged particles in the kinematic range I77I < 2.5 and pj > 100 MeV, 

• at least one charged particle in the kinematic range \t]\ < 2.5 and pj > 2.5 GeV. 

The correction procedures as well as methods used to extract the systematic uncertainties are identical 
to the three phase-space regions presented in the main part of the paper. The first phase-space region 
is chosen to be compared with the other diffraction-reduced phase-space region with six particles above 
500 MeV and allows the study of the interplay between the number of particles and the pj, in particular 
for the study of diffraction models. The second additional phase-space region is chosen so as to be 
less influenced by non-perturbative parts of the non-diffractive modeling and to be useful for predicting 
high-7?T particle rates, for example for trigger studies. 

Table |9] shows the number of selected events and tracks for these two additional phase-space regions 
at both = 0.9 TeV and ^/s = 1 TeV. Figures [15] to [18] show the four kinematic distributions. Table [TOl 
shows the results for the mean track multiplicity at central eta (obtained as the average between -0.2 < 
77 < 0.2). Figure [19] shows the mean track multiplicity at central rapidity for all centre-of-mass energies 
and phase-space regions presented in this paper, along with predictions from pythia6 AMBTl. 



Phas 
rich 


e-Space Region 
min pt 


^fs = 
Events 


0.9 TeV 
Tracks 


Events 


= 7 TeV 
Tracks 


20 
1 


100 MeV 
2.5 GeV 


69,833 
19,016 


1,966,059 
22,233 


4,029,563 
1,715,637 


153,553,344 
2,690,534 



Table 9: Number of events and tracks in the two additional phase-space regions and energies considered 
in this appendix. 



Phase-Space Region 


Energy 
(TeV) 


dnch/di] at ?7 = 
Measured 


rich > 20, ;?T > 100 MeV 


0.9 

7 


6.596 ± 0.025 (stat) ± 0.080 (syst) 
9.077 ± 0.005 (stat) ±0.157 (syst) 


«ch > 1, /^T > 2.5 GeV 


0.9 

7 


0.281 ± 0.006 (stat) ± 0.0005 (syst) 
0.362 ± 0.001 (stat) ± 0.002 (syst) 



Table 10: duchldr] at 77 = for the additional two different phase-space regions considered in this paper 
for = 0.9 TeV and = 7 TeV. 
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Figure 15: Charged-particle multiplicities as a function of the pseudorapidity for events with n^h > 20, 
Pt: > 100 MeV (a,b) and n^h > 1, ;?t > 2.5 GeV (c,d) and \t]\ < 2.5 at ^ 0.9 TeV (a,c) and ^ 
7 TeV (b,d). The dots represent the data and the curves the predictions from dilferent MC models. The 
vertical bars represent the statistical uncertainties, while the shaded areas show statistical and systematic 
uncertainties added in quadrature. The bottom inserts show the ratio of the MC over the data. The values 
of the ratio histograms refer to the bin centroids. 
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n^i, > 20, > 100 MeV, h | < 2.51 
ATLAS = 0.9 TeV 
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Figure 16: Charged-particle multiplicities as a function of the transverse momentum for events with 
rich > 20, pt > 100 MeV (a,b) and rich > I, Pt > 2.5 GeV (c,d) and I77I < 2.5 at = 0.9 TeV (a,c) 
and ^fs - 1 TeV (b,d). The dots represent the data and the curves the predictions from different MC 
models. The vertical bars represent the statistical uncertainties, while the shaded areas show statistical 
and systematic uncertainties added in quadrature. The bottom inserts show the ratio of the MC over the 
data. The values of the ratio histograms refer to the bin centroids. 
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Figure 17: Charged-particle multiplicity distributions for events with «ch ^ 20, /jj > 100 MeV (a,b) and 
«ch > 1, ^-T > 2.5 GeV (c,d) and I77I < 2.5 at ^ 0.9 TeV (a,c) and ^ 7 TeV (b,d). The dots 
represent the data and the curves the predictions from different MC models. The vertical bars represent 
the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties added in 
quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio histograms 
refer to the bin centroids. 
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Figure 18: Average transverse momentum as a function of the number of charged particles in the event 
for events with ?ich > 1, Pt > 2.5 GeV and \r]\ < 2.5 at a/s ^ 0.9 TeV (a) and ^ 7 TeV (b). The dots 
represent the data and the curves the predictions from different MC models. The vertical bars represent 
the statistical uncertainties, while the shaded areas show statistical and systematic uncertainties added in 
quadrature. The bottom inserts show the ratio of the MC over the data. The values of the ratio histograms 
refer to the bin centroids. 
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10^ 10"^ 

\fs [GeV] 

Figure 19: The average charged particle multiplicity per unit of rapidity for ?/ = as a function of 
the centre-of-mass energy. All the measured phase-space regions and energies are shown as triangles 
and compared to predictions from pythia6 AMBTl tune. The phase-space region label is above the 
corresponding curves and points. Combined statistical and systematic uncertainties are approximately 
equal to or smaller than the data points. 
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N. Huseynov^", J. Huston^^, J. Huth^'', G. lacobucci^"^^, G. lakovidis^, M. Ibbotson^^, I. Ibragimov^'^^ 

R. Ichimiya^^, L. Iconomidou-Fayard"^ , J. Idarraga^^^, M. Idzik^^, P. lengo'*, O. Igonkina^*^^, 

Y. Ikegami^^, M. Ikeno^^ Y. Ilchenko^^ D. Iliadis^^"^, D. Imbault^^ M. Imhaeuser^^'^, M. Imorii^^ 

T. Ince^o, J. Inigo-Golfin29, R loannou^ M. Iodice!34a^ q lonescu'^, A. Irles Quiles^^^ K. Ishii^^ 

A. Ishikawa^^, M. Ishino^^, R. Ishmukhametov^^, T. Isobel5^ C. Issever^i^ S. Istinl^^ Y. ltoh^^\ 

A.V. Ivashini28, W. Iwanski3^ H. Iwasaki^^ J.M. Izen'^o, V. Izzoi^^a, b. Jackson^^o, j.N. Jackson^^ 

P Jackson!'^^ M.R. Jaekel^^ V. Jain^\ K. Jakobs'^^ S. Jakobsen^^ J. Jakubek^^^, D.K. Jana^", 

E. Jankowski^^*^, E. Jansen^^, A. Jantsch^^, M. Janus^^, G. Jarlskog'^, L. Jeanty^^, K. Jelen^^, 
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V.A. Ko^otkovl2^ O. Kortner^^ S. Kortner^^ V.V. Kostyukhin^o, M.J. Kotamaki^^, S. Kotov^^ 
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